Shifting donor-acceptor photoluminescence in N-doped ZnO 
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We have grown nitrogen-doped ZnO films grown by two kinds of epitaxial methods on lattice-matched 
ScAlMg04 substrates. We measured the photoluminescence (PL) of the two kinds of ZnO:N layers in the 
donor-acceptor-pair transition region. The analysis of excitation-intensity dependence of the PL peak shift with 
a fluctuation model has proven that our observed growth-technique dependence was explained in terms of the 
inhomogeneity of charged impurity distribution. It was found that the inhomogeneity in the sample prepared 
with the process showing better electrical property was significantly smaller in spite of the similar nitrogen con- 
centration. The activation energy of acceptor has been evaluated to be ~ 170 meV, which is independent of the 
nitrogen concentration. 



PACS numbers: 78.55.Et, 81.15.Fg, 71.35.Cc, 72.15.-v 
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Wide band-gap semiconductors are of intense interest for 
many applications, mainly including optoelectronic devices. 
Many of the wide-gap semiconductors are known to be grown 
more or less under residual n-type even under nominally un- 
doped conditions. The trial of the /5-type doping often results 
in an observation of the radiative recombination from donor- 
acceptor pairs (DAPs) because of the coexistence of these 
two kinds of impurities. The DAP photoluminescence (PL) 
bands have been widely used for more than three decades to 
characterize doped semiconductors. As a result, considerable 
progress has been achieved in description of the DAP band- 
shape and recombination dynamics in a low-doping regime 
(the impurity density being 10^^ cm~^ We are interested 
in the DAP PL for the samples at higher doping levels (in ex- 
cess of lO^'* cm^'^) because such a large amount of nitrogen 
has to be doped to achieve p-type conductivity in the case of 
technologically -promising wide-gap semiconductor ZnO. De- 
spite the tremendous world-wide activities and the successful 
demonstration of homoepitaxial Ught emitters 0], this prob- 
lem has not been solved up to now and remains the major 
bottleneck of ZnO-based optoelectronics. Even when one of 
the best growth techniques for ZnO:N that have ever reported 
[a repeated-temperature modulation (RTM) technique f2\] is 
used, the compensation ratio of p-ZnO is typically as high as 
0.8. 

In our previous work fs*], a theory developed by Thomas, 
Hopfield, and Augustyniak 1 1 ] (THA model) has been 
adopted for the analysis of the PL properties of such samples. 
According to this theory, the PL energy of the DAP recombi- 



nation includes the Coulomb term and is given as: 



hu = hiy^ 



Vei? = [Eg -Ea- Ed) + e^eR, (1) 



*Author to whom correspondences should be addressed; electronic mail: 
[makino @ sci.u-hyogo. ac.jp | 



where the R corresponds to the inter-pair distance between 
donor and acceptor, while the other notations take their con- 
ventional meanings. The activation energy of acceptors (Ea) 
amounted to be as large as sa 260 meV, which is much larger 
than the value (170 meV) obtained spectroscopically by Look 
and coworkers 0], or than the value (« 100 meV) evaluated 
electronically by ourselves f?]. 

A relatively high compensation ratio of ZnO:N films may 
be responsible for a reason of this discrepancy in Ea- Indeed, 
Baume and coworkers reported that their approach based on 
the THA theory failed to obtain good agreement with experi- 
ments for highly-compensated ZnSe:N (5|]. It is thought that 
some of the doped nitrogen atoms act as a donor and the total 
density of donors inside our doped film is rather close to the 
effective density of states for donors (1^x10^* cm^^). In 
that density regime, one cannot neglect effects of overlapping 
of the wave functions of somewhat delocalized donor elec- 
trons. Because the optical transition becomes resembling with 
so-called "free to acceptor" emission, the resulting PL spec- 
trum is then determined by the donor-band position/width and 
the position of the quasi-Fermi level. The corresponding sit- 
uation is schematically shown in Fig. 1. On the other hand, 
in the THA picture, the lineshape of PL spectrum reflects the 
distribution of their inter-pair distances. Moreover, because 
such a high doping level tends to yield a severe inhomogene- 
ity in a statistical distribution of involved charged impurities, 
the importance of the potential fluctuation effects should be 
also accounted for. 

In this work, we show that, by taking the wavefunction 
overlapping and the potential fluctuation effects into consid- 
eration, the earlier-mentioned PL recombination behavior can 
be consistently interpreted. 
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FIG. 1: Contribution of the donor band part for the scheme of the 
spectral band formation. The hatched area corresponds to an output 
spectrum. 



TABLE 1: Sample characteristics of our ZnO:N films. The sample C 
was grown by the RTM technique, while the remainings were grown 
by the conventional method. "HITAB" and "conv." stand for high- 
temperature aimealed buffer and conventional growth technique, re- 
spectively. 



Sample 


Growth temperature °C [N], 10^**cm"^ 


Buffer 


Method 


A 


650 2 




conv. 


B 


550 10 




conv. 


C 


470 7 


HITAB 


RTM 



Our epilayer samples were grown on lattice-matched 
ScAlMg04 (SCAM) substrates by laser-beam-assisted 
molecular-beam epitaxy. For nitrogen doping, a radical 
source was operated. The characteristics of our samples have 
been compiled in Table I. It is now well-known that nitrogen 
cannot be incorporated into ZnO at high growth temperatures, 
at which the high crystallinity growth is ensured for the 
undoped case. To solve this dilemma, we have developed 
the RTM to satisfy both high-crystallinity and high density 
of nitrogen. In other words, we repeated a growth sequence, 
in which ZnO:N layer is deposited at low-temperature, 
followed by rapid ramp to high-temperature, and growth 
of high-quality layer with smaller amount of nitrogen. 
The RTM epitaxy allowed one to obtain p-ZnO when the 
nitrogen concentration was set to [N] w 2 x 10^" cm~^ j^l- 
Secondary-ion mass spectroscopy (SIMS) was used for 
the determinations of the total N concentrations (hereafter 
[N]). The PL of the samples was excited by a continuous- 
wave He-Cd laser (325 nm). Maximum density of excitation 
was w 5 kW/cm^^. 

The evolution of low-temperature (5 K) PL spectra with the 
excitation intensity for the sample A is shown in Fig. 2(a). 
The sample was grown without adopting the RTM nor insert- 
ing the annealed buffer layer. The DAP emission band shifts 
to higher energies with increasing excitation intensity. Such a 
shifting PL has been intuitively explained in terms of the nar- 
rowing of average inter-pair distance of the DAPs at higher 
laser influences. 
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FIG. 2: Excitation-intensity-dependent PL spectra taken at 5 K in 
nitrogen-doped ZnO films of the sample A and C. One-phonon 
replica (DAP-LO) is also observed for them. The emission energy 
of the D'AP shows hardly any shift. 
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FIG. 3: The dependence of peak positions on normalized excitation 
intensity for samples A-C. Symbols are the experimental data and 
solid lines are the calculated results according to Eq. (2). The pa- 
rameters used in the calculation are: Rs = 100 A for sample C; 
iis = 35 A for sample A; and _Ra = 15 A for sample B. Also shown 
in subfigure (b) is the expanded scale presentation for the D'^AP. 



The symbols (closed circles and diamonds) of Fig. 3(a) 
show the excitation-intensity-dependent PL peak energies of 
the DAP obtained for our ZnO:N (sample A and B). The PL 
peaks are excitation-intensity-sensitive in the case of sam- 
ples |3| A and B. The observed shift per two-decades of /exc 
are wl8, and w45 meV. One could learn from these results 
that the heavier doping lead to the more excitation-intensity- 
sensitive blueshifted PL. 

Changing the doping and growth methods (sample C) while 
keeping the dopant concentration approximately unchanged 
from that of 'A' yielded an excitation-intensity-dependent PL 
significantly different from the above-mentioned results. The 
traces of Fig. 2(b) and closed squares of Fig. 3(a) correspond 
to its data. The near-band-edge spectrum [Fig. 2(b)] is dom- 
inated by the two kinds of DAP emissions at 3.218 eV and 
3.235 eV. Because the energy separation of these two DAP 
peaks (ss 17 meV) is close to the difference in the shallow- 
and deep-donor activation energies |7], these peaks are as- 
signed to be shallow- and deep-donor-related DAPs (D^AP 
and D'^AP). It is apparent that the D'^AP peak shifts in 10 meV 
with the two-decade increase in excitation intensity, while 
the D'^AP is insensitive to the excitation intensity. Despite 
the fact of similar density of nitrogen dopant ([N], see Ta- 
ble I) with that of sample A, the peak energies for the control 
samples grown by the conventional method are relatively red- 
shifted and sensitive to the excitation intensity change. The 
growth-technique dependence of the shift is interesting to be 
discussed. 

Recently, Kuskovsky et al. have developed a model f§\ to 
quantitatively interpret the PL of the doped materials in terms 
of ionic-charge distribution inhomogeneity and the wavefunc- 
tion overlapping effects. As was schematically shown in 
Fig. 1, this theory assumes change in the quasi-Fermi level 
of the donor dopants (/id) under excitation intensity is related 
to the observed PL peak shift and calculates self-consistently 



the quasi-Fermi levels as a function of the magnitude of the 
fluctuations, (related to the screening radius Rg) and the exci- 
tation intensity. The calculation showed, in the "fluctuation" 
model, if the donor impurity band is sufficiently naiTow, the 
change of the quasi-Fermi level has almost no influence on 
the resulting PL peak energy. It should be noted that the ef- 
fect of the moving acceptor quasi-Fermi level is much weaker 
because its effective density of states is typically significantly 
greater. 

Quantitatively, the luminescence intensity Isig) per unit 
energy at an energy i? as a function of excitation intensity g 
can be written as follows. We here use E = hv — Eq + 
Ea + Ed as an emission energy instead of hi^. It is important 
to realize that while hv is always positive, E can be either 
positive or negative, depending on the magnitude and a sign 
of the fluctuation energy term. 

Using freezing temperature Tg, Eq. (2) gives the spectrum 
of the zero-phonon line of DAP PL in the presence of fluctua- 
tion as a function of excitation intensity. 
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with. 



u{u + 1 — exp[— u]) 



P{u) = l[l + nC u /e] 



(3) 



Here, we used the normalized excitation intensity of g — 
—2y/TrgiJ,ave- The value Wq is related to the transition prob- 
ability of DAP and erf[a;] is the error function. It must be 
also explained about the other notations used in Eq. (2): 



V 



/(ikBTgeRs); 



E = eRsE/e^; ^ = 2R,/Rb 
n = ttINa - Nd\R%; and fi = iiD^/eRs/2e'^Tg. 

Though there are three running parameters {Rs, Tg, and 
Na — Nd) included in this model, we keep them fixed except 
Rs (which is intimately related to the strength of fluctuations 
and to the width of the donor band) for the present calculation 
to clear up the role of Rs- In all cases, we used 3.182 eV for 
{Eg - Ea - E'dX which corresponds to = 3.428 eV, 
Ea = 171 meV PI, and Ed = 75 meV IIH. It should 
be noted that the energy position of the bound exciton emis- 
sion remained almost unchanged from the value of undoped 
films, probably suggesting the absence of the band-gap nar- 
rowing induced by the nitrogen doping. The solid lines in 
Fig. 3(a) show the calculated results based on Eq. (2) as a 
function of the normalized excitation intensity whose range is 
10^^ < g < 10^"^ for three different screening radii {Rs)- 
It should be noted that the progression of Rg from sample A 
to sample C is reasonable, i.e., smaller Rs for severer red- 
shifted emission. The emission energy {hv) becomes nearly 
independent of the excitation intensity if Rs becomes greater 
than 150 A. The calculated result for the D'^AP of Fig. 3(b) 
using Ea — 171 meV and Ed — 60 meV (corresponding to 



(2) 



4 



the shallower donor) shows small shift, though the averaged 
peak energy is in reasonably good agreement with the exper- 
imental value. It is likely that Rs for D^AP is greater than 
that of D'^AP, as large as ~150 A, probably because of the 
difference in their densities. Although the quantitative agree- 
ment between theory and experiment is beyond our scope of 
work, all the experimental values for samples A-C could ride 
on positions close to the corresponding theoretical curves as 
shown in Fig. 3(a) in the pumping range from g w 10"^ to 
g w 10^"*. Value of g has been evaluated by integrating a 
theoretical spectrum and by some averaged energy on quasi- 
Fermi level 1 8]. We must point out here that the effects of fluc- 
tuations {Rs) play relatively important role for the observed 
difference in the excitation-intensity-dependent shift as well 
as the redshifted PL energy in the samples grown by two dif- 
ferent methods. Using Ejj — 60 meV for the cases of sam- 
ples A and B evidently requires the reduced Rs for the better 
matching with the experiment. But such a change results in 
an excitation-intensity dependence, much more drastic than 
the experiment. 

The same conclusion in terms of inhomogeneity could be 
drawn by comparing the PL widths for three samples studied 
in this work. The present analysis offers alternative and com- 
plementary estimate of the fluctuations than the analysis based 
on the width evaluation of the DAP bands because superim- 
posing D^AP contribution tends to prevent from the precise 
evaluation of the broadening in D'^AP. 

By applying the fluctuation model, we found that the ef- 
fect of inhomogeneity can well explain the growth-technique 
dependence of the donor-acceptor pair (DAP) PL properties. 
The analysis allows the more precise evaluation of Ea and 
also tells us that the improved growth technique can decrease 
the potential fluctuation contained in the sample. The nitro- 
gen atoms are likely to be regularly incorporated by faithfully 
following their energetics under this condition because the 
growth of ZnO:N can initiates on a nearly-step-free smooth 
surface, which presumably results in a narrower distribu- 
tion of impurities. To check the reproducibility, we have 
also grown several RTM samples in different temperature se- 
quences, the dopant density of which are comprable with that 
of the sample C. Their PL peak energies are hardly dependent 
on the exciting laser flux. 

In the pioneering work on this theory, experimental data 
taken for three ZnSe:N samples have been also present in 
which their peak shifts are comparable with respect to each 



other f?], which in turn lead to: 24 A < i?s < 35 A. It 
is worth while mentioning that our present experimental work 
could verify that the theory is more widely applicable for large 
range of fluctuations (15 k < Rs < 100 A). It might be 
interesting to compare [N] evaluated by the SIMS with the 
screening radius Rg. These two quantities link with respect to 
each other in the following formula: 

In the samples A and B, the deduced density from the above 
equation is comparable with or larger than [N], while the den- 
sity for the sample C is smaller than it (probably corresponds 
to the density of D'* or D^). 

It can be proposed that the low-temperature simple PL tech- 
nique combined with this fluctuation model provides a quick 
tool for the assessment of the charged impurity distribution 
randomness and for the evaluation of Ea, compared to the se- 
lective PL spectroscopy or scanning capacitance microscopy, 
both requiring more complicated instrumentation . 

In summary, we have grown ZnO samples doped with 
nitrogen and studied their PL spectra. The RTM-grown 
ZnO:N sample showed both D^AP and D'^AP bands with 
zero-phonon energies of 3.218 and 3.235 eV, respectively. We 
also discussed the PL shifting behavior and compared them 
with the controls grown by conventional method 1 3]. We have 
explained how to extract the degree of dopant distribution in- 
homogeneity from the results of optical characterization. The 
activation energy of acceptor has been evaluated with the use 
of proper modeling to be « 170 meV, which is independent 
of the nitrogen concentration. 
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